Conditional SIC-POVMs 

Denes Petz0, Laszlo Rupperil^ and Andras Szanto [f] 



Department of Mathematical Analysis, 
Budapest University of Technology and Economics, 
Egry Jozsef u. 1., Budapest, 1111 Hungary 

^—t . Abstract 
O; 

. In this paper we examine a generalization of the symmetric informationally 

I complete POVMs. SIC-POVMs are the optimal measm'ements for full quantum 

■ tomography, but if some parameters of the density matrix are known, then the 
^ ■ optimal SIC POVM should be orthogonal to a subspace. This gives the concept 

I of the conditional SIC-POVM. The existence is not known in general, but we give 
■ a result in the special cases when the diagonal is known of the density matrix. 
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>! 1 Introduction 

^■ 
[<. ■ 

i/-^ • The motivation for a positive operator valued measure (POVM) is in the quantum infor- 

^ ■ mation theory. The outcome statistics of a quantum measurement are described by (one 

O , or more) POVMs. A sequence of measurements on copies of a system in an unknown 

£^ [ state will reveal the state. This process is called quantum state tomography |10j . 



A POVM is a set {Ei : 1 < i < k} oi positive operators such that Yli Ei = I. A 
^ ■ quantum density matrix p can be informed by the probability distribution {Tr pEi : 1 < 

^ ■ i < k}. A density p G M„(C) has — 1 real parameters. To cover all parameters k > 

should hold for the POVM. We can take projections -Pj, 1 < « < n^, such that 

VP, = nJ, TrP,P, = ^- (t^j), E, = -P, 
^-^ n + 1 n 

1=1 

and this is called symmetric informationally complete POVM (SIC POVM) by Zauner 
and it is rather popular now [H El IH El [H [H] • Zauner shoved the existence for n < 5 
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and there has been more mathemtical and numerical arguments [6l [16] . The existence 
of a SIC POVM is not known for every dimension. Another terminology for this is tight 
equiangular frame. We may also consider less than rt^ projections with similar properties. 

A SIC POVM {Ei : 1 < z < n^} of an n-level system is optimal for several arguments. 
For example, the SIC POVM was optimal in our paper [12] where the minimization of 
the determinant of the average covariance matrix was studied. Actually, this kind of 
optimization is too complicated and a different argument was in [15], minimization of 
the square of the Hilbert-Schmidt distance of the estimation and the true density. In 
the present paper the minimization of the square of the Hilbert-Schmidt distance will be 
used. 

In this paper the subject is the state estimation again, but a part of the ri^ — 1 param- 
eters is supposed to be known and we want to estimate only the unknown parameters. 
A POVM {Ei : 1 < z < A;} is good when k < rt^ and — k parameters are known. It is 
obvious that the optimal POVM depends on the known parameters and we use the ex- 
pression of conditional SIC POVM. This seems to be a new subject, the existence of such 
conditional SIC POVM can be a fundamental question in different quantum tomography 
problems. The description of the conditional SIC POVM is the main result in Section 1, 
however, the existence is not at all clear. The formalism is in a finite dimensional Hilbert 
space l-i and a state means a density matrix in Bil-L). The known parameters determine 
a traceless part B C B[l-i) and the operators of the conditional SIC POVM are orthog- 
onal to B. In Section 2 a particular situation is studied, we assume that the diagonal 
entries of the state space are given. A mathematical subject called planar difference set 
in projective geometry is used there. 

2 The optimality of conditional SIC-POVMs 

We examine the case of M„(C). Let us suppose that ai is an orthonormal basis of 
self-adjoint matrices, i.e. 

(^i = (^h = kj^ ^, j e {0, 1,2, . . .n^ - 1}. 

We fix (To = -^In- (The elements of this basis are often called generalized Pauli matrices.) 

A quantum state p satisfies the conditions Tr p = 1 and p > 0. It can be written in 
the form 

n2-l 
i=0 

where 9o = A necessary condition for the coefficients can be obtained: 

E^' = Trp^<l- (1) 

i=l 
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We decompose M„(C) to three orthogonal subspaces: 



MniC) = A®B®C, (2) 

where A := {AJ„ : A G C} is one dimensional. Denote the orthogonal projections to the 
subspaces A, B, C by A, B, C. A density matrix p G M„(C) has the form 

p=— + Bp + Cp. 

n 

Assume that Bp is the known traceless part of p and Cp is the unknown traceless part 
of p. We use the notation p^< = p — 'Bp. The aim of the state estimation is to cover p*. If 
the dimension of B is m, then the dimension of C is — m — 1. For the state estimation 
we have to use a POVM with at least N = — m elements. To get a unique solution 
we will use POVM with exactly N elements: {Fi, F2, . . . , F^}. For obtaining optimal 
POVM, we will use similar arguments to [11] which was a straightforward extension of 
the idea appeared in [T5] . 

If {Qi : 1 < i < N} are self-adjoint matrices satisfying the following equation 

1 ^ 
p* = -/ + ^ OiCi = ^PiQi, Pi = TrpFi, 

then {Qi : 1 < z < A} is a dual frame of {F^ : I < i < N}. Then the state 
reconstruction formula can be written as 

N 

1=1 

We define the distance as 

N 

\\p* - P*ll2 = Tr (p* - p^f = {p{i) -p{i)){p{j) -p{j)){Qi,Qj) 
and its expectation value is 

N 

N N N 

i=l i=l j=l 

N 

= Y,V{^{Q^,Q^)-^^{p*f■ 

i=l 

We concentrate on the first term which is 

N 

J2{TTF,p){Q„Qi) (3) 

i=l 
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and we take the integral with respect to the Haar measure on the unitaries l]{n). 
Note first that 

UPU* d^iiU) 



'U{n) 

is the same constant c for any projection of rank 1. If XlILi ~ then 



nc 



It n 

and we have c = Inl'n- Therefore for A = XliLi ^i-^i have 

r " / 

/ UAU* d^iiU) = V A,c = — Tr A 

and apphcation to the integral of (E]) gives 

/ Ti Fi{UpU*)dfi{U) = -TiFi. 
J n 

So we get the following quantity for the error of the state estimation: 

1 ^ 

T := E {\\Up*U* - Up*U*\\^^ dp{U) = - ^(Tr F,)(Q„ Q.) - Tr (p*)^ 

1=1 

This is to be minimized. Since the second part is constant, our task is to minimize the 
first part: 

TV 



5^(TrF,) (Q„Q,) (4) 

We define the superoperator: 

F = J2\F.,){FMTTFr'. 



i=l 



N 



i=l 



It will have rank A^, so if < the inverse of F does not exists, but we can use its 
pseudo-inverse F~, so F~|(Tj) = 0, if cTi G B. Ri is the canonical dual frame of Fi, if 



1^.) =F-|P.), 



where Pi = (TrFi)-^^^. 



Lemma 1 For a fixed Fi, ^ is minimal if Qi = Ri, i.e. if we use the canonical dual 
frame. 
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Proof. Let us use the notation Wi = Qi — Ri. Then 

N N N 

Y,'^^Fi\Ri){Wi\ = 5^TrF,|i?,)(g,| -^TrF,|i?,)(i?i 

i=l i=l i=l 

N N 



J2 Tr F,F- (Q.I - 5^ Tr F,F- |P,) (P,|F- 



i=l i=l 

N N 



F- ^ Tr P, I P,) (g, I - F- ( 5^ Tr I P,) (P, I ) F- 

1=1 i=l 

F n- F FF = F n- F n = 0, (5) 



where 11 = A + C, and we use that from 

TV 

IP 



i=l 

follows 

N 

n = Y,m{F,\. 

i=l 

So we have 

N N N 



Y,T^TF,{Q„Qi) = 5^TrP,(iy,,iy,) + 5^TrPi(iy„P,) 

1=1 i=l i=l 

N N 

+ ^ Tr P, (P„ ly,) + ^ Tr Fi (P, , Ri) 

i=l 1=1 
N N 

= 5];TrP,(iy„W,) + ^TrP,(P„P,) 

1=1 i=l 

N 

> ^TrP,(Pi,Pi). 



i=l 



□ 



We know the optimal dual frame for a fixed POVM Pj, and the following lemma 
provides a property for the optimal POVM: 

Lemma 2 The quantity in ^ is minimal if 
Proof. From ([5]) we have 

N 



5^(TrP,)|P(z))(P(z)|=F-n = F- 



i=l 



so we have the equation: 

N 



^(TrF,)W),i?(^))=Tr(F-). 



i=l 

Let 1^1 > > • • • > t'n2 be the eigenvalues of F. Since the rank of F is A^, we 
1^1 = for i > N. We want to minimize 

^ 1 

Ti.(F-)=5:-. 

i=l 

It is easy to check that A is an eigenfunction of F with ui = 1 eigenvalue: 

N N N 

F\I) = 5^(TrF,)|P«)(P«,/) = ^(TrF,)|P(.)) = = 1^) 

i=l i=l i=l 

and we have the following condition: 

N N N 

Y^u, = TrF = ^(P„Pi)TrF, < J^TrP^ = Tr J = n. 

i=l i=l i=l 

Combining these conditions we get that the measurement is optimal if z/2 = 1^3 = . 

- jn- 

Now we can obtain that the optimal POVM is a conditional SIC-POVM: 
Theorem 1 // 

then 

N 

Vp, = -J, TrP.P, = — — — (z^j), TrafcP, = K G P). 
^-^ n n(N — 1 

1=1 ^ ' 

Proof. Let us use notation Aj = TrPj, then has the form: 

5^A.|P.)(P.|=A + ^C. 



i=l 



Then we have to the following equation: 

N _ . 



i=l 
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with Q := Pk -d- I. 

From {Pi\Q) = Tt PiPk — d the left hand side of ([7]) becomes 

N 

>^^{Q\P^){P^\Q) = h{l - df + P^^k " d) 
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We can compute the right hand side as well: 

A{Pk - dl) = APk -dl = A{Pk - I/n) + I/n - dl = I{l/n - d), 

{Q\A\Q) = {1/n - d)Ti {Pk - dl) = n{l/n - df 
When Pk = Y.f=oCiai, then 



C\Q)=Y^ c,a., {Q\C\Q) = 



cl 



So ([7]) becomes 

Afc(l -dy + J2 ^^(Tr P.Pfc - df = n{l/n - df + ^ cl (8) 

From ([1]) we have 

J2c^<l-cl = l- 1/n. (9) 

This implies 

Tl — 1 

Afc(l - df < nil/n - df + ^^^^(1 - ^/n), 

which is true for every value of d, so 

. n{l/n-df + ^^il-l/n) 

Afc < mm — 

d i^-df 

By differentiating we can obtain that the right hand side is minimal if: 

N -n 

d 



n{N-l) 
and then we get 

n 

Since A^ = n, we have Ai = A2 = . . . = A^f = n/N. 

From that follows that there is an equality in ([9]) too, so we have 

^ c,^ = 1 - Cg d = 0, if ai e B Tt aiPk = 0, if (Ti e B. 



On the other hand from ([8]) we have 
So it imphes 



N — n 

Ti P,Pk = — r if t^k. 



□ 



One has to be careful about this result though, since we only consider the case of 
linear state reconstruction, as it was stated in [15]. Finding the optimal statistic in a 
more general setting requires complicated nonlinear optimalization. 

Now we look at some examples related to the previous theorem and we take different 
values. 

Example 1 If we do not have any information a priori about the state (m = 0, = n^), 
then 

^ n+1 ^ 

so the optimal POVM is the well-known SIC-POVM (if it exists [S]). 

Example 2 If we know the off-diagonal elements of the state, and we want to estimate 
the diagonal entries (m = n"^ — n,N = n), then from Theorem [T] it follows that the 
optimal POVM has the properties 

n 

Tr PiPj = (i^j), ^Pi = I, and Pi is diagonal. 

i=l 

So the diagonal matrix units form an optimal POVM. □ 



Example 3 If we know the diagonal elements of the state, and we want to estimate the 
off-diagonal entries (m = n— l,N = n'^— n + 1), then from Theorem [1] it follows that 
the optimal POVM has the properties 

^-^ n 

1=1 

and Pi has a constant diagonal. More about this case is in the next section. □ 
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3 Existence of some conditional SIC-POVMs 



Theorem [T] tells that conditional SIC-POVMs are the optimal measurements if they exist, 
but it was not written anything about the existence of such POVMs. The existence of 
SIC-POVMs for arbitrary dimension is not known and they are a special case of the 
conditional SIC-POVMs. We can not expect to give a full description of SIC-POVMs, 
but this section contains a particular example. There are seqveral equiangular frames 
with less than projections |3], but it is not clear, what parameters are spanned by 
their complementary part, ie. what the known parameters are. Intuition suggests that 
the case when the known part corresponds to a subalgebra of the full matrixalgebra is 
especially interesting. 

Suppose we know the diagonal elements of a n-dimensional density matrix. We 
want to construct the related conditional SIC-POVM, that is subnormalized projec- 
tions Pi forming a symmetric POVM and complementary to the diagonal projections 
Ei = |ej)(ej| G M„(C) (1 < i < n). These projections form a maximal abelian subalge- 
bra. Easy dimension counting shows, that we want to construct N = — n + 1 such 
projections. 

So {|ej) : 1 < z < ra} is an orthonormal basis in the space. We set 

n 

V " ■ 1 

and a diagonal unitary 

f/ = Diag(g"^g"^g"^...g""), 

where the integer numbers < < — 1 are differents. Another unitary T permutes 
the eigenvectors of U : 



|ej+i) if 1 < z < n — 1, 
|ei) if z = n. 



Note, that T|0) = T*|0) = We have 

|(f/V,e,)P = |(0,(t/*)S-)l' = k~'"^fl(0>e,)r = |(0,e,)p 

= |(0,T^--V)p = |((T*rV,ei)P = |(0,ei)r 

and the projections Pk '■= (p) (U^ (f)\ are complementary to the diagonal projections: 

Tr|t/V)(t/V|(|ei)(ei| -J/n) = 0. 

It is easy to check that 

N N N 

k=l ^ fc=l ^ k=l ^ 
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so we obtain 



N 



^-^ n 

k=l 

and the sum is multiple of I. 

We need to choose the numbers ai, 02, . . . , such that 



TtP.P, = \{U'^\W4>)\' = - 



m=l 



is constant when i ^ j. From the formulas 



5^TrP,P, = (iV-l)lt + l, 5^TrP,P, = Tr (p.J^P,) 
j 3 \ j / 



N 
n 



we obtain t = n — 1. 

Next we use a terminology from the paper [7]. The set G := {0, 1, . . . , — 1} is an 
additive group modulo A^. The subset D := {aj : 1 < z < n} is a difference set with 
parameters (A^, n, A) when the set of differences — aj contains every nonzero element 
of G exactly A times. When this holds, then we have 



1=1 



AT-l 



m{ai-aj) 



+ Ag^ = n - A, 



s=l 



where g is from f jTOj) . Here A = 1. If the appropriate difference set exists, then there 
exists a conditional SIC-POVM. Similar constructions of tight equiangular frames related 
to difference sets are examined in detail in [8]. 

The existence of difference sets with parameters {N, n, 1) is a known problem, named 
the prime power conjecture [171 U\j and we get the following result: 

Theorem 2 There exists a conditional SIC-POVM with respect to the diagonal part of 
a density matrix if n — 1 is a prime power. Then N = — n + 1 and the projection Pi 
(1 < i < N) have the properties 



N 



i=l 



N 



n 



n — 1 



^ 3) 



A few examples about M = {a^ ■ k} is written here: 

n = 2, M = {0,1}, n = 3, M = {0,1,3}, 
n = 4, M = {0,1,3,9}, n = 5, M = {0, 1, 4, 14, 16}. 
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